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Abstract. Several newly developing hybrid imaging methods (e.g., those combining 
electrical impedance or optical imaging with acoustics) enable one to obtain some 
auxiliary interior information (usually some combination of the electrical conductivity 
and the current) about the parameters of the tissues. This information, in turn, 
happens to stabilize the exponentially unstable and thus low resolution optical and 
electrical impedance tomography. 

Various known instances of this effect have been studied individually. We show 
that there is a simple general technique (covering all known cases) that shows what 
kind of interior data stabilizes the reconstruction, and why. Namely, we show when the 
linearized problem becomes elliptic pseudo-differential one, and thus stable. Stability 
here is meant as the problem being Fredholm, so the local uniqueness is not shown 
and probably does not hold in such generality. 
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1. Introduction 

The fast developing in the recent years hybrid imaging methods refer to a range 
of techniques in medical imaging in which different modalities are used in concert 
to benefit from the strengths of each while mitigating their individual weaknesses 
[llIZlEllSniETllSlllHOllSl] . For example, ultrasound tomography provides high resolution, 
while not necessarily providing high contrast. On the other hand. Electrical Impedance 
Tomography (EIT) and Optical Tomography (OT) can provide high contrast but are 
typically plagued by high instability and thus poor resolution [T5l[T8l|22]. Acousto- 
Electric Tomography (AET), also called Ultrasound Modulated Electrical Impedance 
Tomography (UMEIT) [l[5l[26l[30l[32l[33ll57], uses focused (physically or synthetically) 
ultrasound waves to perturb the conductivity inside the object of interest, and these 
perturbations can be measured using the EIT techniques. A similar approach is 
used in UMOT (Ultrasound Modulated Optical Tomography), which combines optical 
tomography with acoustics [3l[l3lll0lill|5l]. The MREIT and CDI/CDII techniques 
use various combinations of EIT and MRI [3714391 B2l [53] . What is common for 
these (physically rather different) techniques is that after some manipulations (see 
[lEHHlinHllESlEalETHaniliaES] for details), the values through the interior of the 
object of a function (often of the form cr{x) |Vu(x)|^) can be obtained. Here a is 
the conductivity to be determined, and u is the corresponding electric potential. The 
problem then becomes to determine a from this interior data. 

Surveys of the recent results on hybrid methods can be found, for instance, in 

It has been observed that a common feature of these methods is that they 
provide significantly better stability and resolution than more conventional EIT and 
OT techniques. The opinion has been expressed by several experts that some meta- 
statement should exist that claims that "appropriate" interior information stabilizes 
the exponentially unstable problems such as EIT or Our aim is to provide a 

version (in fact, several versions of various generality) of such a statement. Doing this, 
we can address the stability of different hybrid methods with internal data in a unified 
way. 

Our considerations are local, i.e. in a neighborhood of a known smooth background. 
We thus use linearization. If one could prove that the Frechet derivative of the 
corresponding non-linear mapping exists, has zero kernel, and is a (semi-) Fredholm 
operator in some scale of (say, Sobolev) Banach spaces, by standard functional analysis 
arguments this would have implied local uniqueness and stability of the non-linear 
problem. In this text, we provide very general statements on (semi-) Fredholm property 
of such Frechet derivatives (see, e.g., [27] - [29| 155] for basics on Fredholm and semi- 
Fredholm operators). This is done by reducing them to elliptic pseudo-differential 

X The main difficulty in EIT and OT is that a signal by the time it reaches the boundary has largely 
"forgotten" where it originated from. Thus the idea is that having some quantity F(a{x), u{x), Vu(a;)) 
attached to an arbitrary interior location x should compensate for this "lack of memory." 
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operators. This achieves necessary stabihty estimates modulo the possible finite- 
dimensional kernel. The authors doubt that the uniqueness claim (i.e., absence of 
non-trivial kernel) can be made in such wide generality in which we obtain ellipticity 
and thus Fredholm property. Thus, the absence of kernel should be dealt with on case- 
by-case basis (as it has been done before). However, existence of Fredholm property 
shows that the interior data does have stabilizing effect on the initially exponentially 
unstable problem. 

Let us describe the structure of the article. In Section [2], we describe the framework 
for later sections, and we state basic definitions and lemmas. In Section [3] we investigate 
inverse conductivity problems with different types of additional interior data. We first 
consider the data of the form a{x)\'Vu{x) |^, which is known to arise in a variety of hybrid 
problems |i[6HHl[IlHlll26l[30l[33l|3M9l|4a[55 For p G (0, 1) we show that ellipticity, 
and thus stability arises with a single set of the interior data. For p > 1, two such 
sets are needed. Finally, we look at a rather general type of interior data of the form 
F{a{x), u{x), Vn(x)) and obtain sufficient conditions under which one gets the problem 
stabilized. In Section H] we treat a problem coming from the so-called Quantitative 
Photo-Acoustic Tomography (QPAT) p^V[T2l[23] - [25l[30l[33] . where the equation becomes 
of the diffusion type (with an absorption term) rather than just a divergence type 
equation as in the previous Section. 

The main results of the paper are contained in Theorems 13.21 13.61 13.81 and 14.11 
Sections [5] and [6] contain the proofs of some technical lemmas used and final remarks 
and conclusions correspondingly. 

2. Preliminaries 

Let f2 be a bounded open region in with smooth boundary, and let Q' be an open 
region compactly contained within Q. We will also need an intermediate domain Q": 

n' (E n" m n. 

We will frequently need to use two cutoff functions, xi ^ C^{Q) that is equal to 1 in a 
neighborhood of fi" and X2 ^ that is equal to 1 in a neighborhood of Q'. 

In Section [3] we will address various inverse conductivity problems, in which the 
goal is to determine the unknown log-conductivity a in the elliptic problem [§| 

L^u = -V ■ (e'Vw) = (1) 

in Q from some boundary data. For instance, in the EIT (Electrical Impedance 
Tomography) the data might be the whole Dirichlet-to-Neumann map on dQ [T8 | [T9 | [22 | 

EEHIH]. 

§ It will be more convenient for us to work with log-conductivity, which we denote a, rather than with 
the true conductivity exp a. 
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We thus will have to work frequently with the Dirichlet boundary value problem 
-V ■ (e'Vw) =0 
u\dn = J- 

It will be important for us that the solution of ([2]) depends on log-conductivity a. To 
emphasize this dependence, we will sometimes write u{a) for the solution to ([2]). 

Analogously to [21], we define the affine space of admissible log-conductivities as 

= I ^\n\n' = 0} . (3) 

Functions in can be considered as defined on by extending them by zero. 

We assume that a E L^di^) and / E H^/^{dn) in so uE H\n) (see. e.g. [H]). 
In Section H] we will study the following more general problem: 

L^^u := -V ■ (e^Vu) + e^u = 

u\dn = J- 

The coefficients a and 7 are the log-diffusion and log-attenuation coefficients, 
respectively. 

Lemma 2.1 The map from [a, 7) to u, defined by Q), is Frechet differentiate as a 
mapping from L°°{n) x to H\n) at any point (ao, 70) E x 

This fact is well known, but we supply in Section [5] its proof. Infinite differentiability of 
the coefficients (cxo, 70) is in fact an overkill assumption here, but we will use (and thus 
prove) the lemma only in the smooth case. 

In many hybrid imaging methods (see, e.g. |i|l5|[7tl26 | [33 |[37] - [39|H2l[53] ). internal 
information of the form 

F((t(x),m(x), Vm(x)), (5) 

where u is the solution of ([2]) corresponding to some specific boundary data /, can be 
derived from the measured data. Thus, the next goal is to recover the log-conductivity 
function <t{x) from the knowledge of F{a{x),u{x), Vm(x)) for all x E Q (which explains 
the name "interior data"). Since, for a fixed boundary Dirichlet data /, the functions u 
and Vu{x) are determined by a, we can consider F as a non-linear operator acting on 
a: F -.a^ F{a). 

The pseudo-differential technique that we use makes the statements simple and 
their proofs rather transparent. We make use of basic facts about pseudodifferential 
operators on (see, for instance, or [l5l Ch. 7]). We will use the standard symbol 
classes ^'"(M") comprising smooth functions a{x,C,) that satisfy for any multi-indices 
a, (3 and for sufficiently large |^| the estimates 

d:d^^a{x,o\<c^p{i + \^\r-^^^ (6) 

with some constants Ca/3- (Such symbol classes are often denoted S'™o(I^"); but we will 
omit the subscripts as we will not be considering more general symbol classes.) The 
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corresponding classes of pseudo-differential operators, which we denote by 0P5''"(]R"), 
are given by 

a{x, D)u = J I a(x, Oe'^^"y>^uiy) dy (7) 
where a{x,^) G S'™(]R") (with the standard regularization of this expression, see e.g. 

Hang). 

If a symbol a{x,0 e 5'"(M") satisfies 

a{x,^) = amix,^) +r{x,^), (8) 

where r(x, ^) G S'"'^^(]R"), we will call am{x,^) the principal symboloi a(x,^) G S'™(]R"). 
The principal symbol is determined modulo 5''"~^(M'^). If, for some i? > 0, the principal 
symbol am{x,S,) satisfies the estimate 

\a^{x,i)\>C\irioi\i\>R, (9) 

then the symbol a{x,^) is called elliptic, and the corresponding operator a(a;, D) is called 
elliptic as well. 

We will also need some facts about (square) matrix pseudo-differential operators. 
Let Ai^j{x,^) for i,j = l,...,p be a matrix of classical symbols of pseudo-differential 
operators. Suppose that there exist two p-tuples {si,...,Sp} and {ti,...,tp} of real 
numbers such that Aij G S''^*"'"*-' (R'^). Let also A^j G 5'**"'"*^ be their principal symbols. 
The system A = {Aij} is said to be elliptic in the Douglis-Nirenberg (DN) sens^Ul (see [1] 
or [m Ch. 7], [52l Ch. 9]) if the determinant det {Al^^x.C)) does not vanish for |^| > i?, 
for a suitable R. 



3. Stability in Inverse Conductivity Problems with Internal Data 

Here, we address functionals F of the more specific form 

F((t) = e"|Vn(a)|^'. (10) 

Several hybrid imaging methods provide this kind of internal data (see, e.g., [UE]- 
[8l[nHia[26l[33l[3M9lll2l[53]). For example, it arises with p = 2 in AET (Acousto- 
Electric Tomography, also called sometimes UMEIT (Ultrasound Modulated Electrical 
Impedance Tomography), or Impediography) . 

When < p < 2 and a G L^(f2), one concludes that F{a){x) belongs to L^(f2). 
Let us consider F as a non-linear mapping from L^(f2) to L}{VL). As such, it is Frechet 
differentiable at any smooth log- conductivity ctq, as long as the corresponding solution 
u{ao) has a gradient that is bounded below by a positive constant. This is a direct 
consequence of the following lemma: 

II Sometimes it is also ealled Agmon-Douglis-Nirenberg elliptieity. An equivalent, although differently 
formulated, notion was also introduced by L. Volevich (e.g., [52l Ch. 9]). 
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Lemma 3.1 Let F{y, z, w) be a function of three variables that is smooth when y, 2; G M, 
and w > 0. Assume that F satisfies the bound 

\Fiy,z,w)\<Ciy)iz' + w'), (11) 

where C{y) depends continuously on y. Then the mapping a =^ F{a, u, |Vm|) is Frechet 
differentiable at the smooth background log- conductivity ctq, as a mapping from 
to L^ip). 

The proof is provided in Section |5l 

Since the Frechet derivative exists, it can be found by a formal hnearization. 
Consider a small perturbation of a^: 

u = uo + ev + o[e), 

where uq = u^ctq). A simple substitution, as in [21], shows that v G Hq{Q) solves the 
boundary value problem 

L^,v = V ■ (pe-°Vno) ^^^^ 
v\dn = 0. 

One notices that the dependence of f on p is linear. We will indicate it as v{p). 

The Frechet derivative dF of F, which we will denote by dF, is a linear bounded 
operator from L^(f2) to L^[Q). Applying the chain rule and ( l3Ti) to e'^|Vn(cr)|P, one 
finds dF as 

dF{p) = pe-"|Vnor + pe-° ^;^;y . (14) 



|Viio 



We introduce a cutoff version of dF, which extends to a pseudo-differential operator on 
M". Let xi be a smooth cutoff function supported in Q which is identically equal to 1 
on a neighborhood of Q". We define an operator A mapping L°°(f2) to L^(R") by 

A{p) = xidF{xip)- (15) 

Because of the presence of the cutoff by xi before applying dF, the operator A has 
a natural extension to L°°(M"). 

In order to show that A is a pseudodifferential operator, we analyze equation ( 13T]) . 
The expression V ■ (xipe"^" Vmq), as a differential operator acting on p, has principal 
symbol ie'^°Xi^ ' Vuq. This operator, when acting on functions that vanish outside of 
Q', does not depend on the choice of Xi- 

The principal symbol of the elliptic differential operator L„g is e^^^l^p. Hence, L^^ 
has a pseudo-differential parametrix P G OPS'~^(M") with principal symbol (e'^'^ |^p)~^. 
This means that L^rgP = I + S, where / is the identity operator on M" and S is 
a smoothing operator on M", and analogously for PLo-g. Let us define the following 
function: 

u;:=P(V-(xipe"«Vuo)). 
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Then we have 

L^^{v -w)=V- (pe'^^Vno) - L,,P{V ■ (xipe'^" Vuo)) 

= V ■ (pe'^" Vmo) -{I + S){V- (xipe'^" Vmo)) 

= V ■ ((1 - xi)pe""Vno) + 5(V ■ (xipe'^" Vuo)). (16) 

The expression in (fT6l) is a smooth function, so La^^{v — w) E C^^Q). By elhptic 
regularity, v = w mod C°°(fi). Because of this equivalence, the mapping p ^ v is 
a pseudo-differential operator modulo infinitely smoothing operators on Vt. All other 
operations in equation (fT4|) are simply multiplication operators, so we see that after 
multiplying by Xi? ^ is a pseudodifferential operator on W^. 

Let Aq{x,^) denote the principal symbol of A, so that A{x,^) = Ao{x,^)+R^i{x,^) 
where -R_i is a symbol of order —1. The symbol is easily derived. The principal 
symbol of a composition of operators is the product of the individual principal symbols. 
Applying this to the composition of the operators 



P ^ V ■ (xipe""Vno) 

shows that the mapping p k-)- v{p) has a principal symbol given by 

(^e-Xi^-Vno) (e-|en-^ = ^^^^|p^. 
From equation (1141) we then find that Aq is given by 



M^,0 = ^i^'l^^or + PX?e-y^|^^^. (17) 



Let 6 denote the angle between ^ and Vuo. Then the principal symbol Aq becomes 

Ao{x,0 = xle'-'lVnoini-pcos^O) (18) 

on Q'. Since we are interested in determining under what conditions A is elliptic, we 
are therefore motivated to consider the case when p < 1 separately from p > 1. 



3.1. The p <1 case 

Theorem 3.2 If p < 1, then 

(i) A{x,D), as defined above, is a pseudo- differential operator of order zero, which is 

elliptic in a neighborhood of Q" ; 
(a) dF , as an operator acting in L'^{Q'), is Fredholm; 

(Hi) Let K be the kernel of dF as an operator on L'^{Q'), and let R be its range (K is 
finite- dimensional and R is of a finite co- dimension). Then dF, considered as an 
operator from L'^{Q')/K onto R, is a topological isomorphism, i.e. there exists a 
constant C such that for p G L^{Q'), 

< IIpII <C\\dF{p)\\ (19) 
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Remark 3.3 

When we consider dF as an operator from L^(f2') to itself, we are really considering 
the operator T o dF, where T : L^(M"') — )■ is the restriction operator. The 

same goes for A. 

((i) Sometimes we need to consider T as acting into L^(]R"), in which case the action 
of T is simply multiplication by the characteristic function of Q' . The operator T 
also maps the corresponding Sobolev spaces: 

T : i7^(M") ^ H'{Q'). (20) 

Proof. Since Aq{x,C,) = Xi^^^" I Vuo|^(l —pcos^O) and p < 1, A is elliptic on a 
neighborhood of Q". This proves the first claim of the theorem. 

Let X2 be a smooth cutoff function supported in the fl" that is equal to 1 on a 
neighborhood of Q'. We define a symbol Q{x,^) G S'°(M"') by setting 

Then, letting I be the identity operator, Q{x, D)A{x, D) — X2{x)I is a pseudodifferential 
operator of order —1 on MJ^. 

We claim that TQ is a left inverse for TA modulo compact operators on 
Indeed, 

TQTA -I = {TQA - I) + TQ{T - I) A. (22) 

Since the pseudo-differential operator Q{x, D)A{x, D)—X2{x)I is of order —1 on M", the 
function {TQA — I)f belongs to H^iVL') for any / G L^(f2'). In addition, the function 
(T — I)Af is equal to zero on so by the micro local property of pseudodifferential 
operators, TQ(T — I)Af is a smooth function on Q'. Thus the righthand side of equation 
(!22|) is a bounded operator from L'^{Q') into H^{Q'). The imbedding operator of H^{Q') 
into L^(f2') is compact, proving the claim. 
We also observe that 

A{x, D)Q{x, D) = Q{x, D)A{x, D) (23) 

modulo operators of order —1, so TAQ — I maps L^(f2') — )■ H^{Q'). Therefore TA is also 
right-invertible modulo compact operators. This means that TA is a Fredholm operator 
on L^(fi'). Then T o dF is a Fredholm operator on L'^{Q') as well, because dF = A as 
operators on L^(f2'). This proves the second claim of the theorem. 

The 3rd claim of the theorem is an immediate consequence of the 2nd one. □ 

3.2. Thep>l case 

lip > 1, the formula (fTSi) shows that there are directions ^ at each point x in which 
the principal symbol Aq{x,$,) = Xie"^" I Vuo|^(l —pcos^O) vanishes. In order to make 
the problem elliptic, we need to assume availability of two measurements. Namely, for 
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two different boundary conditions /i, /2 in ([2]), let u^'^^ for i = 1 , 2 be the corresponding 
solutions: 

L^u^'^ = 
u^^lan = U 

Assume that we are given knowledge, for all x G f2, of the following sets of internal data: 
F22 = e''" IVu^^^l" (25) 

Such functionals have been extracted from the measured data in hybrid imaging 
methods, see for example [8l[2T|[33]. 

Our (vector) internal measurement function will be now 

/ ^ii(^) \ 

F{a):= F22(a) . (26) 
V ^12 (^) / 

We again consider small perturbations of a smooth background log-conductivity 
(To as in f lT2|) . and let Uq^ be the corresponding solutions of with ctq as the log- 
conductivity: 

r L^.ul!^ = 

(27) 



uyidn = Ji- 

We need to assume in addition that the gradients of Mq^"* and Mq^^ are nowhere parallel: 
4'^ (28) 
This is known to be possible p] in 2D under an appropriate choice of /i and f2- However, 



as shown in [3S], in 3-D, it is not always possible to choose boundary conditions such that 
( 128|) is satisfied. Condition ( 128|1 means in particular that the gradients are nonvanishing. 

Lemma 3.4 

(i) As in the previous sub-section, functionals Fu and F22 are Frechet differentiable 
with respect to a as mappings from to . 

(a) When p < 2, the mapping F12 is Frechet differentiable if u'q^ and u^q^ satisfy 

■ Vu'-o^ > a (29) 

for some a > 0. // this condition fails at some points, it is still true that the 
functional (f)Fi2 is Frechet differentiable if (p is a smooth cutoff function and [29\) 
is satisfied on its support. 

(Hi) When p = 2, mapping F12 is Frechet differentiable without condition ^29) . 
Remark 3.5 
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• We will handle the case 1 < p < 2 and note that the case p = 2 follows by the same 
argument, with the part concerning the smoothness of F12 omitted. 

• For the case p = 1, the functional F12 is simply not needed, though the following 
arguments still apply. Thus the full range 1 < p <2 will be covered. 

The proof of this lemma can be found in Section [5l 

Given a vector ^, let 61 be the angle between ^ and 'Vu'^\ 62 be the angle between 
^ and Vug^'', and 6 be the angle between Vmq^'' and Vu^q \ As in the case when p < 1, 
we define cutoff versions of dF^ by setting Aij = xidFijXi- The principal symbols of 
All and A22 near Q', calculated in the same manner as before, are 



An(x,0 =X?e'^lv4'^r(l-pcos2^i), 
A22{x, = xy^lVu'i^ fil-p cos^ 62). 

The principal symbol of A12 at points where Vu^q^ and Vmq^^ are not perpendicular is 
also easily derived from the formula 



dF, 



12 



(2) 



+ e 



p 



p/2 
p/2-1 



■ Vt/'^(p) + Vu'o' ■ Vv^'^ip) ) (30) 



,(2) 



(2) 



,(1), 



with 



V ■ (pe'""VMo) 
/*^|an = 0, 

which readily follows from formal linearization. Hence, 



Ai2{x,0 



p/2 



p/2 



COS' 



p COS 61 COS 02 

I COS 6 1 



(31) 



.(32) 



If, near each point {x,C,) E ^l' x ]R"\0, at least one of these symbols is non- vanishing, 
then the (vector) operator {An A22 ^12)* is overdetermined elliptic. 

We notice that Aii{x,^) vanishes when cos6'i = Near points (x,^) where 

both Aii{x,^) and ^22(3;, ^) vanish, the symbol Ai2{x,^) will have to save the situation. 
Near such a point, i.e. where | cos6'i| = | cos6'2| = a simple trigonometric estimate 
shows that the expression | cos6'| is separated from zero, and thus the symbol Ai2{x,C,) 
is smooth. The non- vanishing of Ai2{x,C,) then boils down to cos 6' 7^ 1 at the point x, 
which is guaranteed by ( l28l) . 

Let us define three symbols homogeneous of order on f2 x M", ipii[x,^), ip22{x,C,), 
and ipi2{x,^), such that 

^(X, = V^ll(x, OAllix, + V^22(X, 0^22(X, + ^12(X, O^^X, (33) 

is a non- vanishing symbol of order near Q' x M"\0, and therefore bounded away from 
zero by the compactness of the cosphere bundle of Q'. The above arguments imply 
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that ipi2 can be taken to be zero near points where condition ( I29p is not satisfied. (For 
p = 1 we do not need the functional F12, so we set iIji2 = in this case.) The operators 
ijjij{x,D) with symbols tpijlx,^) are bounded on L'^{fl). This means that the operator 
norm of \Ef(x, D) can be controlled by the sum of the operator norms of Aij{x, D). We 
thus have 

Theorem 3.6 

For 1 < J5 < 2, the operator 

dF : L^{n') {L^n')}"^ 

is semi-Fredholm with a possible finite dimensional kernel K . 
((i) Letting K be the finite- dimensional kernel of dF , the estimate 

\\p\\L2{n')/K ^ IM-^(p)ll{L2(n')}3 - ^ \\p\\L^n')/K (34) 
holds for p G L^{Q') and some constant C . 

Remark 3.7 ^45 we have already mentioned, such functional arise naturally and have 
been studied previously. When p = 2, a similar local stability estimate was proved 
in I21\ \33 f in the space C^'°'{Q'), and in a global estimate was established in W^'°°{Q). 
A stability estimate for a single functional of the form e'^\Vu\'^ was also established in 
on a part ofQ. In the case p = 1, inversion procedures and reconstructions for a single 
functional were obtained in f37\^3^. 



Proof. As in the proof of Theorem 13.21 let 

Then the operator Q{x,D) lies in 0P5''^(]R"). Arguments identical to the ones in the 
proof of Theorem 13.21 show that T\l/(\E'ii, \E'2, "^u) is left regularizeil^ for the operator 

(Fu\ 

F22 : L\n') ^ {L\n')Y. (36) 

Thus, the operator in (1361) is semi-Fredholm with a finite dimensional kernel (and 
infinite-dimensional co- kernel), which implies all the statements of the theorem. □ 



3.3. More general interior data 

We next consider a single, rather general functional and formulate a sufficient condition 
for the corresponding linearized problem being elliptic (and thus Fredholm). 

Let F{y^ z,w) be a function of three variables satisfying the conditions of Lemma 
13. Ij that is, F smooth when y, 2; G M, and w > 0, and satisfies the bound 

% Operator B is a. left regularizer to operator A, if BA — / is a compact operator [^[55] , 
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\F{y,z,w)\<C{y){z' + w' 



where C{y) depends continuously on y. 

The Frechet derivative can be derived by a formal calculation as before, and is given 

by 



, dF dF dFVuo-Vv 



(37) 



As calculated before, the principal symbol of the operator mapping p to f is 
iC, ■ VmqI^I"^, which is of order —1. The middle term on the right hand side in expression 
f l37p is therefore of lower order than the other two terms on the right hand side and 
does not influence the ellipticity of the overall principal symbol. Also as before, we set 
A = XidF. From the other two terms in ( l37|) we find that the principal symbol of A is 



AJx n = ^ - ^iilZ!^!^ 

°^ '^^ dy dw iVnolleP 



— - — IVmoIcos^^ 
oy aw 



(38) 



near Q' (we omit the cutoff function xi because it is identically equal to 1 there). This 
leads to a sufficient condition for the ellipticity of Aq: 



Theorem 3.8 // 

dF{ao,Uo, IVmoI) 



dy 



>\\Vuo\\ 



dF{ao,uo, |Vmo|) 



dw 



(39) 



pointwise in a neighborhood of Q, then 

(i) y4(x, ^) is elliptic of order on a neighborhood of Q' ; 
(a) dF as a Fredholm operator in L^(f2'); 

(Hi) Letting K be the finite- dimensional kernel of dF , the estimate 
1 



C 



IM^(p)llL2(n') < \\p\\l^u')/k < C\\dF{p)\\ 



(40) 



holds for pe L^{n'). 



Remark 3.9 (i) The functional F{a) = (t|Vm|^ satisfies assumption for p < 1, 
so Theorem \3.8i generalizes Theorem \3.2[ 

(a) // (To = c for some constant c, we can have |Vuo| = 1 by selecting an appropriate 



boundary condition, e.g. f = Xi. 

n. 



Then just means 



dF 
dy 



> 



— I 

dw I 



pointwise on 



Proof. We may assume that xi = 1 on the region where condition ( l39l) is satisfied. 
If (139|) holds, then there exists 6 > such that 



dF 




dF 


dy 




dw 



llV^oll 



L°°{n) 



>6>0 



(41) 
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on ft by the continuity of the terms involved. This imphes that Aq{x,$,) > S, and so A 
is elhptic on the same neighborhood of Q'. By setting 

and composing with the restriction operator T, we obtain an operator T o Q[x,D) 
which is a left and right inverse for TA on L^(f2') (and hence for T o dF as well) modulo 
compact operators. Estimate (HOj) then follows. □ 



4. Stability in Quantitative Photo- Acoustic Tomography 

The standard model for diffusive regime photon propagation in biological tissues is 

u := -V ■ (e'^Vn) + e^n = 

(43) 

u\dn = J 

(see, e.g. [51]). The coefficients cr and 7 are the log-diffusion and log-attenuation 
coefficients, respectively. We will assume in this section that a G L^(f2) and 7 G 

The PAT (Photo Acoustic Tomography) procedure done first, provides one with 
the values inside Q of the function 

F{x) = T{x)e^{x)u{x). (44) 

Here, T{x) is the so-called Griineisen coefficieni0 describing the transfer of 
electromagnetic into acoustic energy, which we assume here to be identically equal to 1. 

This function is the initial data for the Quantitative Photo-Acoustic Tomography 
(QPAT), which strives to reconstruct the coefficients a and 7 from the data ( l44l) . 

We will denote by Fj{x), j = 1, 2, . . . , J the internal data (jHj) that correspond to 
solutions of ( l43l) with different boundary conditions fj. 

For a pair of such measurements (Fi, F2), according to Lemma [2.11 the mapping 
((7,7) =^ F := {Fi, F2) is Frechet different iable at smooth background coefficients 
(<^0)7o)- The derivative can be computed formally as before: 



a = (To + ep 

7 = 7o + 

= 4^') + ev^^^ + o(e) (45) 

where z/, p G Lg°(f2). Substitution into (l43l) shows that v^^^ G H^{^1) solves the boundary 
value problem 

Sih _ ,,^7o„,(i) 

(46) 



+ The Griineisen coefficient is in principle also not known, so one might want to include it as an 
unknown in the reconstruction procedure, e.g. |12j . We are not doing this here. 
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We thus find that the differential of the mapping Fj is 



14 



(47) 



A calculation similar to the one in the previous section shows that the operator 
Lcro^'yo has a parametrix with principal symbol 

1 



(4J 



e^o|^|2 + e7o 

which is equivalent to (e'^°|^p)~^ modulo lower order terms. Hence, according to (|4 
the matrix of the principal symbols of the operator (p, z/) i— >■ xidFxi is given by 



A{x,0 ■=Xi 



151 = 



u. 



(2) 



(1) 



(2) 



Xi 



(49) 



modulo lower order terms. 

We consider here the Douglis-Nirenberg parameters 

S = {Sl, S2) = (-1, -l),t= {h, t2) = (0, 1) 

and attempt to check the DN ellipticity of dF. 

The determinant detA{x,^) is non- vanishing near Q', if 

e ■ (n^^Vn^^) - n^^Vn^^)) ^ . 



(50) 



Thus, ellipticity fails at the vectors ^ orthogonal to the field 
The natural idea is to have more measurements that would provide a basis of vector 
fields and thus preserve ellipticity. This leads to the following result: 

Theorem 4.1 Suppose that one has access to 2n measurements 

(-^1,1; -^1,2), • • • , (-^n,l, -^n,2), 

such that the vector fields 

for k = 1, . . . , n form a basis at each point x in Q. We define the operator F as follows: 

F := (-Fl,l, -Fl,2, • • • , Fn^i, Fn,2) ■ 

Then the operator 

dF : L\n')^H^{n') {H\n')Y 



I 2n 



is semi-Fredholm with a finite dimensional kernel. 

Letting K be the (finite- dimensional) kernel of dF , the estimate 
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1 
C 




< 



H/K 




< C 



_H'i(f7")2" 




(51) 



H/K 



holds for some constant C > 1. 

Here we used the shorthand notation 



Remark 4.2 The assumptions we made on the vector fields Vk, as the reader could see, 
arose naturally. It is interesting to notice that they are the same that were also arising 
in the study of QPAT in [T4\l , in which the CGO (complex geometric optics) solutions 
technique was used. The authors of derive a global estimate that is somewhat similar 
to (EIP; 



\\5a\ 



c\n)' 



c\n) 



< C 



E 

k=l,...,r. 
J=l,2 



|-Ffc,i(0-1, 7l) - Fk,j{(T2, 72) 



(52) 



for I > 2, where 5a = ai — a^-, ^7 = 7i — 72- 

Under an additional convexity assumption on dQ, the authors of [T^ also derive 
such an estimate with only two vector fields, when I >3. 

The reader might also notice that in l[51\) . in comparison with ^5M}, different orders 
of smoothness are used, which allows us to get the two-sided estimate. 

Proof. The 2n x 2 matrix operator X2dFx2 has the principal symbol 



A{x,0 ■=X2 



(1,1) \ 



(1,2) 



V 



u 



{n,l) 



U 



{n,2) 



X2- 



J 



Here the principal symbol is understood in the Doughs- Nirenberg sense with parameters 
s = (—1, —1, . . . , — 1, — 1), t = (0, 1), and X2, as before, is a smooth cutoff function that 
is equal to 1 in a neighborhood of f2' and vanishes outside il". 

Under the assumptions of the theorem, at every (x,^) G fi' x (R"\0) the symbol 
A{x,^) is injective (since at least one of the square 2x2 blocks 



u, 





(fc,2) 



for k = 1, . . . ,n is invertible). Thus, the operator is overdetermined elliptic in DN sense. 

Thus, there exists a left parametrix with the principal 2 x 2n symbol B{x,^) with 
DN parameters s = (0, — 1), t = (1, 1, . . . , 1, 1) (this is a well-known construction, which 
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we indicate at the end of Section E]). In other words, for the corresponding pseudo- 
differential operators dF and B one has Bx2dFx2 = I + K m a, neighborhood of VL' , 
where K is a. smoothing operator. 

Let us agree to extend both p G L'^{Q!) and v G Hl{Q!) as equal to zero outside ^2', 
without changing notations for these extended functions. Then, according to the DN 

parameters, X2dFx2 ^ ^ ^ belongs to (iJQ(f2))^", and the corresponding mapping from 

Hq{Q') to (i^Q (fi))^" is continuous. This proves the right hand side inequality 
in dSU. 

Similarly, B acts continuously from (//^^^^(M"))^" to Ll^^p{W) ^ H^^^piW) . 
Then for the restriction of I ^ | to fi' we can write 



X2Bx2dFx2 \ ]-X2K' 



or 



(/ + X2K) 1^1= X2BdF 



(;) 



This gives us the estimate from above of the following kind: 

+ X2K) iy^J^ \\L^{n')®Hl(si') < C\\dF ^ ^ I ||/fi(n")2-- 

Since the operator X2K is compact m. H = L'^{Q") ^ Hq{Q"), the operator I + X2K is 
Fredholm in H. 

This implies the remaining statements of the theorem. □ 

5. Proofs of some lemmas 

After proving Lemma 12. 1^ we then proceed to give a proof of Lemma 13.11 Finally we 
prove Lemma 13. 4[ making use of Lemma 13.11 

5. 1 . Proof of Lemma \2.1\ 

Proof. First of all, we reduce (j4]) to a problem with homogeneous boundary conditions. 
Let E the operator of harmonic extension from dVt to Vt. Then, replacing the solution 
u with t>o-,^ := u — Ef, we reduce (jl]) to 

La,'^V(j,^ = /(T,7, ^^^^ 

Va,')\dn = 0, 
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where 

■■= V ■ (e'^'VEf) - e^'Ef G H~\Q). (54) 
Now the map (a, 7) ^ Mo-,7 factors as the composition of the following maps: 

(a, 7) ^ {(e^ e^), U,^} ^ L^,, ^ L~l^ ^ L^l^U, + Ef. (55) 



The first map in (155|) . is clearly Frechet differentiable as a mapping from L°°{Q) x 
L°°(f2) to L°°(f2) X L°°(f2) x H~^{^1). Indeed, fa,-y depends linearly and continuously 
on (cr,7). 

The second map {e'^,e'~') Lo-,^ is linear and continuous from L°^{Q) x L°^[Q) to 

L{H^{n),H~\n)). 

The operator -^0-0,70 ^ L{Hq{Q), H~^{Q)) is invertible (see the simplest case of this 
statement in [15], Ch. 5, Proposition 1.1] and general results in [T6l[36]); Lcr,^ is thus 
invertible in a neighborhood of (7o,cro) in L°°(f2) x L'^[Q). 

The mapping 

L{H^{n), H^\n)) L{H~\Q), H^{n)) (56) 

of taking inverse operator is known (e.g., [55]) to be analytic on the domain of invertible 
operators. This implies differentiability of the last two mappings in (l55l) and thus proves 
the lemma. □ 

5. 2. Proof of Lemma \3.1\ 

Proof. The function F{a{x),u{x),\Vu{x)\) lies in L^{Q). Indeed, we have 
|F((T,M, |Vm|)| < C(| l2,oo(n))('"^ + l^""!^)) l^""! both square-integrable 

functions. 

As a result of Lemma 12.11 the map 

(7 (a, u, Vn) (57) 

is Frechet differentiable from L^(fi) L'^di^) x L^{n) x {L^{n)}''. (The middle space 
could be taken to be H^{Q) as before, but we will not need this here.) We claim that 
the map 

L^(fi) X L\Q) X {L^Q)}" L\Q) 
{f,g,h)^F{f,g,\h\) 

is Frechet differentiable at {fo,go,^o) ^ -^S(^) ^ L'^i^) x {L^(fi)}" when /o, go, and 
ho are bounded and smooth with |ho(x)| > m for x & Q and some positive constant m. 
The boundedness of /o, go, and ho implies that for any multiindex a, the function 

D^'Fiy, z, w)\y=f^,i^^)^ z=go[x), w=\ho(x)\ (58) 



is a bounded function on f2. 
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Let (/, (y(,h) be a triple of functions in L^{^) x x {L^(n)}"'. Consider the 

function E = E{x) defined on Q by 



E = F{fo + f,go + 9, |ho + h|) - F(/o, go, |ho|) 

- VF(/o,^?o,|ho|)- 

|ho| 

= Fif,g, |h|) - F(/o,^?o, |ho|) - |^(/o,^?o, |ho|)/ 

- ^(/o,^o, |ho|)^ - ^{fo,9o, |ho|)^r^ • 
oz aw |ho| 

The function E lies in since each individual term does. 

We estimate the L^-norm of E as follows: let 

U = Ixen 



max{|/(x)|, \g{x)\, |h(x)|} > m| . (59) 

On Q\U, we can apply Taylor's Theorem to F to find that 

\E{x)\<C{f{xf + g{xy + \h{x)\') , (60) 

where C depends on an upper bound for the second order partial derivatives in f l58|) . 
As a result, 



[ \E{x)\dx<C I (f{xf + g{xf + \\i{x)f)dx 
Jn\u Jn\u ^ ' 

< C(No\{^) ll/llioc(f^) + ||^||i2(J^) + ||h||^2(^)„ 



On [/, we have 



/ \E\<C I (|(7o + ^7r + |ho + h| 



'u Ju 

+ \F{fo,go, IhoDI + |VF(/o,^?o, |ho|)|(|/| + \g\ + |h|) ) . (61) 



The constant in front depends only on F and | |/o| Iloo(-q), as we need to consider only 
1 1/1 li:,oo(Q) < m, say, which makes the constant C{y) in ( ITT]) bounded. Using the fact 
that /o, go, and ho are bounded functions, we have 



[ \E\<C [ (l + |/| + |^| + |^|2 + |h| + |h|2) 
Ju Ju 

<c((l + ll/llLoo(n))Vol(f/) 

+ (ll^?llL^(n) + l|h||^.(n).)Vol(f/)^/2 



(62) 
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Owing to the inequality 

Vol(f/)<^(||^7||^.(^) + ||h||i.(^)„) (63) 
for ||/||ioc(n) < m, we have 

f \E{x)\ dx<C + \\g\?mn) + WH^n)-) ' (^4) 

This proves the Frechet differentiabihty of the map (/, g, h) ^ F{f, g, |h|) at (/o, go, ho). 
Hence by the chain rule for Frechet derivatives, F{a, u, |Vu|) is Frechet differentiable as 
a function of a at ctq. □ 

5. 3. Proof of Lemma \3.4\ 

Proof. The Frechet differentiability of Fn and F22 follows from Lemma 13. 1[ Only 
Frechet differentiability of F12 needs to be proven. 
As a result of Lemma 12.11 the map 

^0 

is Frechet differentiable from L'^{Q) to L^(r2)^" at (Tq. We claim that the map 



.-I-?., I (65) 



I ^ (t)\vi ■ V2\P^^ (66) 

V2 J 

is Frechet differentiable at a pair of smooth vector fields vi and V2 that satisfy 

\vi\, \v2\ < M,\vi- V2\>a > (67) 

on the support of </> for some M > 1. To show this, let wi, W2 G L^(f2)" and define a 
function E G L^(i7) by 



E = <P(^\{VI + Wi) ■ {V2 + W2)f^ - \Vi ■ V2YI'' 

- 2 1^1 ■ ^'2r^^"^(^'l ■ V2){V2 -Wi+Vi- W2)J . (68) 

In order to estimate the L^(r2)-norm of E, define a set U by 



a 



U=^xeQ\m^x\w,{x)\>—j . (69) 
On Q\U, both wi and W2 are bounded above by a/4M, so we have 
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\{vi + wi) ■ {v2 + W2)\ > a- 


1^1 


■ W2\ - 


\V2 


■ Wi\ - 


\Wi ■ W2 


>a- 




a 




a 




\vi 




\V2 




{ " 




AM 


AM 


\am) 



a 

> - 

- A 



(70) 



It therefore follows from ( 1681) and Taylor's formula applied to the function of 2n variables 
\x ■ yp^^ (which is smooth and has bounded derivatives when the arguments satisfy 
a/A <x,y<M) that 



^ \E{X)\ dX<C + ||W2||i2(f^)n) 



'n\u 

where C depends on m, M, p, n, and a. 
On U, we use the triangle inequality: 



(71) 



\E\ < \vi ■ V2 



|P/2 



Vi ■ W2 + V2 ■ Wi + Wi ■ W2 



Vi ■ V2 

+ ^l"^! ■ ^'2r^^~"^K'2 ■ Wi + Vi- W2\ . 



p/2 



We make use of the inequality 

\\l + z\P/^-l\ <Cp\z\ . 



(72) 



(73) 



to estimate the term in the first line of fl72l) . (For \z\ < 1/2 this inequality follows from 
the fact that the function z i— ?■ |1 + has bounded derivative for \z\ < 1/2, while for 
\z\ > 1/2 it follows easily from the fact that p/2 < 1.) After using (1731) . we obtain 



E{x)\ dx < Cp\vi ■ V2Y^'^ ^1^1 ■ W2^ V2 - \ W\ - W2\ 



u 



+ ■ ^'2^^^ ^\V2 ■ Wi + Vi ■ W2\ . 

By the Cauchy-Schwarz inequality and the bounds on vi, V2, and vi ■ V2, 



(74) 



E{x)\ dx < C I \Wi\ + \W2\ + \Wi ■ W2\ 

u Ju 

1/2 



< C(lkillL2(f/)n Vol(f/)i/^ + ||'u;2|li2(f;)„ Vol(f/)i/2 

+ II'"^i|Il2((7)" Il'"^2||i2([7)n ) • 

The volume of U can be estimated by 



(75) 



Vol(f/) < Vol 

<2('Ji 



KI>j^})h-Voi({ 

2 



|^i'2| > 



a 
AM 



I 1 1 2 II 1 1 2 

F'i|Il2(Q)" + IF2||l2(q)„ 



(76) 



Stabilizing Inverse Problems 



21 



Hence, 

/ rfX < C + ||W2||i2(0)n) (77) 

This proves that the map ( l66l) is Frechet differentiable. Since condition (1671) holds for 
the vector fields Vuq^^ and Vuq^^ on the support of 0, the Frechet differentiability of F12 
at (Jo follows from the chain rule. □ 

5.4- Left parametrix for over- determined DN elliptic operator 

We provide here a sketch of the classical construction of the left parametrix B used in 
the proof of Theorem 14.11 

According to the assumptions of Theorem 14.11 there exists a finite open covering 
[Uj] of the compact subset Vt' x S*""^ in Vt" x (M"' \ {0}), such that for each j, there 
exists a number kj such that A^.i^x,^) is invertible for (a;,^) G Uj. Here S"""^ is 
the unit sphere in the ^-space M". Consider a smooth partition of unity tpjlx,^) on 
Q' X S*""^ subordinated to the covering. We can always assume that it is positively 
homogeneous of order zero with respect to C, outside a neighborhood of the origin and 
smooth on Q" x M". Let us also denote by Pk : C^'" t— ?■ the operator such that 
Pfc(ai, a2, . . . , a2n-i, a2n) = (02^-15 Ci2k)- One sees that = PkA. Then one can check 
that the symbol B{x,^) = J2ji^j{^^0^kj{x,C,)~^Pkj is the one we require. 

6. Final remarks 

(i) As we have mentioned, some of the models discussed in this article have been 
studied previously. Reconstructions from the functionals ( p5|) with p = 2 have been 
performed for instance in [HIEIIISS], with similar stability estimates being obtained 
in [3133], the approach in [8] being global. Global reconstruction from a single 
functional was considered in [6] by solving a Cauchy problem inward from parts of 
dQ for a nonlinear hyperbolic equation. In this case, stability results were obtained 
in two dimensions, and in parts of Q in higher dimensions. 

(ii) The p = 1 case has been studied in [37)439] . In [38], an iterative reconstruction 
procedure was provided, whose effectiveness was demonstrated in numerical 
experiments. Our analysis in terms of pseudo-differential operators implies that 
although the inversion of the corresponding linearized operator from a single 
functional is not an elliptic problem, ellipticity only fails at points (x, ^) where 
Vmo II ^- Since most singularities of p would likely not satisfy that condition, we 
should expect accurate reconstruction of sharp features in almost all the 
numerical experiments demonstrate. 

(iii) While we have shown infinitesimal Fredholm property of the problems with internal 
data, infinitesimal uniqueness has not been shown and we suspect that it does 
not hold under our very general conditions. However, uniqueness should hold 
generically, which we plan to address elsewhere. 
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(iv) Besides absence of an infinitesimal uniqueness result, there is another obstacle for 
obtaining the local uniqueness and stability result for the non-linear problem. 
Namely, Frechet differentiability is proven in worse function spaces than the 
Fredholm property. We also plan to address this discrepancy in a future work. 

(v) We suggest that parametrices constructed in this paper could be used for 
approximate reconstructions and for pre-conditioning iterative methods. 

(vi) In the cases when our ellipticity analysis asks for multiple measurements, this does 
not mean that reconstruction with a smaller number of measurements is impossible. 
On the contrary, such reconstructions have been achieved by solving hyperbolic and 
degenerate elliptic problems in [71[32],[33] . However, such approaches naturally lead 
to correspondingly error propagation and some blurring effects for the parts of the 
wavefront sets where ellipticity is lost [30|[33]. 

(vii) The types of the internal data functionals F considered in this paper do not cover 
all the needs of hybrid methods. For instance, some non-local functionals of a arise 
in UOT [3]. We plan to address those in a subsequent work. 
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